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Abstract 

We discuss how much we can probe the effective number of neutrino species N v 
with cosmic microwave background alone. Using the data of WMAP, ACBAR, CBI 
and BOOMERANG experiments, we obtain a constraint on the effective number 
of neutrino species as 0.96 < N v < 7.94 at 95 % C.L. for a power-law ACDM flat 
universe model. The limit is improved to be 1.39 < N v < 6.38 at 95% C.L. if we 
assume that the baryon density, N v and the helium abundance are related by the big 
bang nucleosynthesis theory. We also provide a forecast for the PLANCK experiment 
using a Markov chain Monte Carlo approach. In addition to constraining N u , we 
investigate how the big bang nucleosynthesis relation affects the estimation for these 
parameters and other cosmological parameters. 



1 Introduction 



Cosmology is now becoming a precision science, and cosmo logical observations can give us a 
lot of information for our understanding of the universe. Moreover, the interplay between 
cosmology and particle physics in various contexts has also been discussed vigorously. 
One of such examples is the effective number of neutrino species N v . Although collider 
experiments such as LEP have measured the number of light active neutrino types to be 
2.92 ±0.06 [1], it is important to cross-check this value because cosmological measurements 
may lead to different value. This could be due to an extra radiation component which 
is predicted by some models of particle physics such as sterile neutrinos (see Ref. [2] and 
references therein), or due to incomplete thermalization of neutrinos in the low-scale re- 
heating universe in which the reheating temperature T re h can be as low as T rc h ~ 0(1) MeV 
and N v is predicted to be less than three [3-6]. If such a non-standard ingredient exists, it 
can affect big bang nucleosynthesis (BBN), cosmic microwave background (CMB), large 
scale structure (LSS) and so on; thus precise cosmological observations can probe these 
scenarios through the effective number of neutrino species. 

Constraints on N v have been investigated in the literature using the information of 
CMB and LSS, sometimes with priors on the Hubble constant, cosmic age and Helium 
abundance [7-30]. Although CMB in general can constrain various quantities severely, 
since the effects of N v on CMB are degenerate with some cosmological parameters, the 
studies so far have combined CMB data with some other observations such as LSS to obtain 
a sensible constraint on N u . However, when one uses the data from LSS, constraints can 
become different depending on how one treats non-linear correction/bias on small scales 
for the matter power spectrum [23]. Furthermore, different LSS data seem to give different 
constraints on N v [16,17,20,21,23]. Regarding the prior on the Hubble constant H , as is 
summarized in Ref. [20], it can yield some constraints on N v when combined with CMB 
data (without LSS data) [8-10,12,13], but they depend on the H prior adopted. One 
may consider that we can use the usually assumed prior on the Hubble constant based on 
the result by Freedman et al. H = 72 ± 8 [31], but another group reported a somewhat 
lower value as H = 62.3 ± 5.2 [32]. If the lower value for H is adopted as the prior, a 
resulting constraint on N v would be different. Having these considerations in mind, it is 
desirable to investigate a constraint on N v without these kind of uncertainties. 

In this paper, we study a constraint on N u from CMB experiments alone. By making 
the analysis of CMB data alone, we can avoid such subtleties as the galaxy-bias/non- 
linear corrections and the value for the prior on the Hubble constant. However, as is 
mentioned above, the effects of N v are strongly degenerate in CMB with other cosmological 
parameters such as energy density of matter, the Hubble constant, and the scalar spectral 
index, and, in fact, we could not obtain a meaningful bound only with WMAP3 [20,23]. 
Recent WMAP5 alone analysis gives a better constraint but it still cannot give an upper 
bound [27,28]. As we will discuss later, the degeneracy is significant up to about the 
2nd/3rd peak of the CMB power spectrum where the observation of WMAP has precisely 
measured. To break this degeneracy to some extent, it would be helpful to have the 
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information at higher multipoles where signals unique to relativistic neutrinos are expected 
to appear [33] . Recently, the data from ACBAR which probes CMB at higher multipoles 
than those of WMAP has been updated [34]. By using this data in addition to other small 
scale observations such as BOOMERANG and CBI, we can obtain a relatively severe 
constraint on N u which is comparable to that have been obtained previously with LSS 
data. 

The organization of this paper is as follows. In the next section, we start with the dis- 
cussion how N u affects the CMB power spectrum, which helps to understand our results for 
the constraint on N v . In Section [3l we study the current constraint on N u using observa- 
tions of CMB alone. We use the data from WMAP5, the recent ACBAR, BOOMERANG 
and CBI. Furthermore, we forecast the constraint from the future Planck experiment. In 
the final section, we summarize our results and discuss its implications for some models 
of particle physics/the early universe. 



2 Effects of N v on CMB 



The effective number of neutrino species N v represents the energy density stored in rela- 
tivistic components as 



Prad = P 7 + Pv + Px 



7 /4\ 4/3 

1+ 8 n» v 



p 7 (i) 



where p 7 , p v and p x are energy densities of photons, three species of massless active 
neutrinos and some possible extra radiation components, respectively. In this paper, we 
assume that neutrinos are massless and have no chemical potential. For the case with 
the standard three neutrino flavors without an extra relativistic component, the effective 
number of neutrino is N v = 3.046 where some corrections from the incomplete decoupling 
due to a slight interaction of neutrinos with electrons/positrons and finite temperature 
QED effect to the electromagnetic plasma are taken into account [35]. Any deviation of 
N v from this value implies that there exists an extra relativistic component and/or some 
non-standard thermal history takes place such as the low reheating temperature scenario. 

To illustrate the effects of N v on CMB, we plot CMB power spectra for several values of 
N u in Fig. [TJ Other cosmological parameters are assumed as the mean values of WMAP5 
alone analysis for a power-law flat ACDM model. As seen from the figure, as we increase 
the value of N u , the height of the 1st peak is strongly enhanced and the positions of 
acoustic peaks are shifted to higher multipoles. Also, the amplitude on small scales (higher 
multipoles) is suppressed. Below, we discuss where these changes are coming from. 

One of the main effects of N v comes from the change of the epoch of the radiation- 
matter equality. By increasing (decreasing) the value of N u , the radiation- matter equality 
occurs later (earlier). Thus the increase (decrease) of N v gives an almost the same effect of 
the decrease (increase) of energy density of matter. One of noticeable features is that the 
height of the first acoustic peak is enhanced by increasing the value of A,,. This is due to 
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Figure 1: CMB power spectra for the cases with N v = 1 (blue dotted line), 3 (red solid 
line) and 5 (green dashed line). Other cosmological parameters are taken as the mean 
value from WMAP5 alone analysis for a power-law flat ACDM model. 



the early integrated Sachs- Wolfe (ISW) effect in which fluctuations of the corresponding 
scale, having crossed the sound horizon in the radiation-dominated epoch are boosted by 
the decay of the gravitational potential. Thus a larger amount of relativistic species drives 
the first peak higher. Another effect is the shift of the position of acoustic peaks due to the 
change of the radiation-matter equality through the change of N„. The position of acoustic 
peaks is well captured by the so-called acoustic scale 6 a which is inversely proportional to 
the peak position and written as 

Va = —, r \L) 

ro(z xec ) 

where rg(z rec ) and r s (z rec ) are the comoving angular diameter distance to the last scattering 
surface and the sound horizon at the recombination epoch z iec , respectively. Although 
^e(^rec) almost remains the same for different values of iV,,, r s (z rec ) becomes smaller when 
N u is increased. Thus the positions of acoustic peaks are shifted to higher multipoles 
(smaller scales) by increasing the value of N v . Furthermore, since the position of the n-th 
peak can be roughly written as l n ~ rnr /6 a, separations of the peaks become also greater 
for larger N„. 

Another important effect is free-streaming of ultrarelativistic neutrinos [33]. The per- 
turbation of ultrarelativistic neutrino propagates with the speed of light, which is faster 
than the sound speed of acoustic oscillations of photon-baryon fluid. The coupled photon- 
baryon component behaves to oscillate like a compressional fluid; on the other hand, 
ultrarelativistic neutrinos free-stream to erase their fluctuations. These two components 
are coupled via gravity; thus the fluctuations of photons can also be affected by the free- 
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streaming of neutrinos, which results in the damping of the amplitude and the shift of 
the acoustic oscillations. The effects are significant on small scales where fluctuations of a 
given scale enter the horizon while the energy density of ultrarelativistic neutrinos takes 
a significant portion of that of the universe. 

Although the effects of the standard cosmological parameters on the heights and po- 
sitions of acoustic peaks are well known, here we discuss them in some phenomenological 
way including the effects of N u . For this purpose, we calculated the response of the heights 
and positions of the acoustic peaks to the change of the cosmological parameters up to 
the 5th peak around the fiducial values. As a fiducial parameter set, we take those of the 
mean value from WMAP5-alone analysis for a power-law flat ACDM model. The shifts of 
the positions of acoustic peaks are 

Au b Au m An s AY V Ah AN U , . 

Ah = 15.58 - 26.99 + 36.01 + 0.94—^ - 44.58— + 15.53—^, 3 

u b u rn n s Y p h N v 

Al 2 = 62.57^ - 74.90^ + 14.69^ + 2.81^ - 108.60^ + 47.73^, (4) 

u b u m n s Y p h N v 

Al 3 = 74.23^ - 143.47^ + 9.88^i + 4.64^ - 152.39^ + 81.82^, (5) 

to b uj m n s Yp h N u 

Al 4 = 110.84^ - 181.89^ + 7.29^ + 6.77^ - 220.21— + 112.76^, (6) 

uj h uj m n s Y p h N u 

Al 5 = 136.88^ - 237.18^ + 6.20^ + 7.79^ - 276.02^ + 145.34^, (7) 

u b uj m n s Yp h N v 

where U means the position of the i-th peak. u b , uj m , n s , Y p and h are energy densities 
of baryon and matter, the scalar spectral index, the primordial helium abundance and 
the normalized Hubble constant. n s is defined at the wave number k = 0.05 Mpc -1 . The 
positive derivatives of the peak positions with respect to N v demonstrate the decrease in 
r s (zrcc) due to the increase in N u . 

The responses of the heights of acoustic peaks to the change of various cosmological 
parameters are 

^1 = 0.429^ - 0.632^ - 0.947^ - 0.0065^ + 0.141^, (8) 

C h Lu b iu m n s Yp N v 

AC,,, Aiu b Aum An, AY„ AN U , s 
— ^ = -0.211 - 0.579 - 0.034 - 0.035 v - + 0.083 -, (9) 

Ci 2 u b u m n s Yp N v 

ACu Au b Au m An s AY V AN U , . 

— - = 0.026 - 0.136 + 0.276 - 0.071—^ - 0.080—^, 10 

C h u b u m n s Y p N v 

= -0.044^ - 0.229^ + 0.587^ - 0.125^ - 0.108^, (11) 

Ci 4 to b to m n s Y p N v 

AC lr Auo b Auj m An, AY n AN, y , . 
^ = 0.149 - 0.006 + 0.776 - 0.172 v - - 0.216 -. 12 

C h u b lo„ l n s Y p N v 
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where Ci = 1(1 + l)Cf T /2ir. Since h gives only a negligible change to the height of the 
peaks, we omit it. By increasing the value of N v , the height of the 1st peak is enhanced due 
to the early ISW effect, as previously mentioned. For the third and the higher peaks, the 
heights are damped more by increasing N u , which is inferred from the negative coefficients. 
This is due to the effect of free streaming of neutrinos [33] . 

In addition, for later convenience, we also show the derivatives of the peak heights 
relative to the first peak height following Refs. [36]. Here, Hi = Ci/Ci for % = 2-5. They 
are useful quantities when we interpret degeneracies since the dependence on the overall 
amplitude is cancelled out. 

AH 2 = -0.291 — + 0.023— + 0.396— - 0.013^ - 0.026—, (13) 

io b uj m n s Y p N v 

AH 3 = -0.177^ + 0.206^ + 0.514^ - 0.028^ - 0.098^, (14) 

u b u) m n s Y p N u 

AH, = -0.102^ + 0.082^ + 0.317^ - 0.025^ - 0.054^, (15) 
u b u m n s Yp N v 

Auo b Auj m An s AK AN U , . 

AH 5 = -0.040 + 0.084 + 0.236 - 0.023—^ - 0.052—^. 16 

u b u m n s Yp N v 



3 Constraint on N u from observations of CMB 

In this section, we present our result for the constraint on N v from CMB alone. First, 
we give some details of our analysis. We use the CMB data of WMAP5 [27,28,37-39], 
BOOMERANG [40-42], CBI [43] and ACBAR [34]. We performed a Markov chain Monte 
Cairo (MCMC) analysis to obtain constraints on cosmological parameters using cosmomc 
code [44] with some modifications which are described in the following. We explore a 
9 dimensional parameter space which consists of cu b , u c , r, 9 S , n s , A s , Asz, Y p and N v . 
Here, u c is the energy density of dark matter, r is the optical depth of reionization, 6 S 
is the acoustic peak scale [45], A s is the amplitude of primordial fluctuations and Asz is 
the amplitude of thermal Sunyaev-Zel'dovich (SZ) effect which is normalized to the Cf z 
template from Ref. [46]. 

As far as CMB is concerned, these cosmological parameters can be considered to be 
independent. However, when we take into account the BBN theory, Y p is determined once 
u b and N v are given. In this case, we should relate these parameters to each other and 
sample an 8 dimensional parameter space. We shall refer to this relation among Y p , uj b and 
N u as the BBN relation. For this purpose, we calculate Y p as a function of u b and A^ using 
the Kawano BBN code [47] with some update in the nuclear reaction network part based on 
Ref. [48]. Such relation is considered in the CMB analyses in Refs. [26,49-52]. In passing, 
we would like to make a comment on the fitting formula for the BBN calculation presented 
in Ref. [53] which has been used in the authors' previous work [50]. We do not adopt the 
formula here since, as we will see later, our MCMC chains for the constraints from the 
present data sets sometimes go to the region beyond the range over which their fitting 
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Figure 2: ID posterior distributions of N u . The red solid line uses WMAP5 alone (with 
Y p = 0.24 fixed) and the other lines use WMAP5+ACBAR+BOOMERANG+CBI with 
different assumptions on Y p . The green dashed line fixes it to be Y p — 0.24, the blue dotted 
line uses the BBN relation to fix Y p from and N v , and the magenta dot-dashed line 
treats Yp as a free parameter. For the analysis with WMAP5 alone, we assumed the prior 
on the cosmic age as 10 Gyr < t < 20 Gyr. 



formula is valid, < N v < 6. For the Planck forecast, where the chains are contained 
in that region, we obtain the same results with the Kawano code and the fitting formula 
of Ref. [53]. Furthermore, we also consider the case with fixing the helium abundance 
to Y p — 0.24 since, in most analyses, the primordial helium abundance is fixed to this 
value. Finally, it should be noted that this BBN relation is not necessarily realized in 
some cases. We can think of more exotic scenarios in which the BBN theory cannot relate 
those parameters. For exmaple, Ub and Y p may vary between BBN and CMB epochs [54] 
or an increase in N u may take place [55]. 

Now, we present our results in order. In Fig. [2j the posterior ID distributions for N v are 
shown for the analysis with WMAP5 alone and WMAP5+ACBAR+BOOMERANG+CBI 
(CMB all). The former is shown by a red solid line. For the latter case, the results for 
different assumptions on Y p are depicted: Y p being fixed as Y p = 0.24 (green dashed line), 
Y p determined from the BBN relation (blue dotted line) and Y p being treated as a free 
parameter (magenta dot-dashed line). Corresponding constraints on N u are summarized 
in Table[TJ We also show 2D contours of 68% and 95% C.L. in the planes of N u v.s. several 
other cosmological parameters in Fig. [3j Table [2] summarizes the derived constraints on 
these parameters. 
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Mean 


68% t 
68% 1 


95%t 
95% 1 


WMAP5 


5.65 


7.88 


9.96 


(Y p = 0.24: fixed) 


3.02 


1.92 


CMB all 


4.24 


5.47 


7.94 


(Y p . free) 


2.03 


0.96 


CMB all 


3.71 


4.80 


6.38 


(Y p : BBN relation) 


2.27 


1.39 


CMB all 


3.89 


4.89 


6.84 


(Y p = 0.24: fixed) 


2.19 


1.28 



Table 1: The mean values, 68% and 95% limits of N v for several current CMB data sets 
and assumptions of Y p . 



As seen from Fig. the likelihood for N v from WMAP alone has irregular shape, far 
from Gaussian. It has the maximum around N v ~ 5, declines slowly as N u increases and 
drop abruptly at N v ~ 9. The abrupt cut can be traced to the prior on the cosmic age 
to which is implicitly assumed as 10 Gyr < to < 20 Gyr in the analysis. In particular, the 
lower limit to > 10 Gyr makes the cut (see Fig [3]). We can regard this prior to be very 
conservative on the observational ground since it is far looser than the astrophysical lower 
bound of the cosmic age e.g. t > 11.2 Gyr (95% C.L.) from the age estimates of globular 
clusters [56]. Moreover, we should include such prior from the practical reason. As can 
be seen by the relatively slow decline of the likelihood for 5 < N v < 9 or the elongated 
contours in Fig [31 the degeneracy is so severe that we cannot produce MCMC chains 
which are well converged within a reasonable time. Although we can formally calculate 
a constraint using this data as shown in Table [H since the likelihood is so irregular, we 
would conclude that it is not meaningful to constrain N v from WMAP5 alone. 

However, it may be instructive to understand how the degeneracies arise in the WMAP- 
alone analysis. As clearly shown in Fig. [31 N v most notably degenerates with uj c and Ho (or 
9 S ). There are also some degeneracies with n s and A s but not as severe as u c or H . These 
degeneracies are understood as follows. First, to produce the same amount of the early 
ISW effect, u c has to be increased as N v increases. It roughly scales as (u b + u c ) oc N v to 
make the matter-radiation equality same. At the same time, under the flatness assumption, 
fl m has to be preserved in order to have the same distance to the last scattering surface. 
Then, since Q m = (ui, + uj c )/h 2 , h has to increase for larger N v . The slight enhancement in 
n s and A s can be attributed to their effects to cancel the suppression around the diffusion 
damping scales due to the increase in N v . A more quantitative argument based on the 
derivatives presented in the previous section may be useful. The degeneracy as regards the 
same matter-radiation equality is given by setting AC^ = to be Au m /uj m ~ 0.2 AN U /N U . 
This is equivalent to Auj c /uj c ~ 0.3 AN U /N U , which roughly gives the slope in the N v -u c 
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Figure 3: The 68% and 95% allowed regions in the plane of N v 
parameters when WMAP5 alone is used with Y p 



v.s. other 
- 0.24 (red solid line), 
WMAP5+ACBAR+BOOMERANG+CBI are used with Y p = 0.24 (green dashed line), Y p 
being fixed by the BBN relation (blue dotted line) and Y p being treated as a free parameter 
(orange and yellow shaded region). 



plane in Fig Using this relation with All = shows the N„—h degeneracy. From 
Eq. ([3]), we obtain Ah/h ~ 0.2 AN U /N U , which appear in the N u -h contour in Fig[3J The 
N u -n s degeneracy is given by further requiring AH 2 = 0. Plugging Auj m /uj m and Ah/h 
in Eq. (fl3|) yields An s /n s ~ 0.05 AN V /N V . This 5% increase in the best fit value of n s for 
AN V = 3 is consistent with the N v -n s contour in Fig[3l Although WMAP has measured 
the CMB power spectrum very precisely, since it is just up to around the 2nd peak, the 
effects of N v are absorbed in the changes of u c , h, n s and A s and we cannot constrain ./V^. 

For a visual illustration of the degeneracy, in Fig. HJ we show CMB power spectra for 
the case with N u = 1 (blue dotted line), 3 (red solid line) and 5 (green dashed line) with 
other cosmological parameters being chosen such that they give the degenerate spectra up 
to the 2nd/3rd peak. We can see that these curves coincide up to the 2nd peak but begin 
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Figure 4: An illustration of the degeneracy of N„ with other cosmological parameters. 
Here the value of the effective number of neutrino are assumed as N v — 1 (blue dotted 
line), 3 (red solid line) and 5 (green dashed line) and other cosmological parameters are 
chosen such that CMB spectra becomes the same as that with the fiducial parameters. 



to separate around the 3rd or the higher peaks. On small scales, the change caused by 
N v cannot be fully canceled just by tuning other parameters. In particular, n s affects the 
spectrum in the whole scales; thus, even if the spectra is almost degenerate up to the 2nd 
peak by tuning the value of n s , it cannot cancel the damping on small scales. 

When we include the data at higher multipoles such as ACBAR, BOOMERANG and 
CBI, the likelihood functions now have a well-behaved peak close to Gaussian and we 
can obtain meaningful constraints. The bound on N u is 0.96 < N„ < 7.94 at 95 % 
C.L. when Y p is treated as a free parameter. At higher multipoles, the free streaming 
of neutrinos damps the spectrum, which cannot be compensated by above-mentioned 
parameters. Hence the degeneracy can be removed to some extent. That is the reason 
why we can have a severer constraint on N u by including the data on small scales. In fact, 
Y p also suppresses the amplitude on small scales via diffusion damping (see Eqs. (TTTT) and 
(fl2l)): thus the constraint on N v slightly changes for different treatments of Y p but the 
differences are very small as seen in Fig. [2] and Table [TJ Since current CMB observations 
on small scales are not so precise, it does not make much difference how we treat Yp. 
Imposing the BBN relation tightens the constraint to 1.39 < N v < 6.38 at 95 % C.L., but 
it is not so different from the Yp-free case. Also, the limit does not differ much even if 
we set Y p = 0.24. Similarly, the estimates for the other cosmological parameters are not 
affected by the assumption on Y p as shown in Table [2j 

Up to now, we have assumed no prior on N u . However, if we consider an extra radiation 
component such as sterile neutrinos and so on, the effective number of neutrino species just 
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parameters 


CMB all 

(Y p . free) 


CMB all 
(Y p : BBN relation) 


CMB all 
(Y p = 0.24) 


WMAP5 
(Y p = 0.24) 


/ i 1 

e s 

T 

n s 

^(10*%) 

Y P 
N u 


u - uzzy -0. 00054 

n 1 OO + 0.018 

u.ioz_ 037 

1.03901HS 
0.088±g;ga 
0.977 + [|-^ 

o ln4 +0.0067 
°- lu ^-0 .0061 

n 9?n +0 - 105 

U.ZZU_ 085 

4.24+^? 

i - U ^-0.34 


0229i+°-°° uw 

U.U^ai_ .00059 

n 1 24+ 0017 
U.l^^_ 027 

1 H/1 1 Q+0.0039 
l.U^lO_ 0051 

o-o88±g : a? 

n O7Ci+0.024 

— U.UZO 

o QO7+0.066 
O.UJ 1 _ 061 

n 256+ - 015 

U.ZdU_ 016 

79+1.07 
°- ' z -1.45 

1 n7+ - 93 

L - u ' -0.33 


02291+° °° 58 

u.uzzyi_ 00053 

n i 97+0.017 

U - 1Z ' -0.029 

1-040210™ 
0.0881S 

3.098+°^ 
3.8911;?° 

i - ud -0.34 


n n997^+ u - ouoti0 

U.U^/J_ .00062 

n 1 e;Q+0.036 
U.±OO_ 047 

1.03341°;°°- 
0.088+°;°!* 
0.989+°°?° 

q i 90+0. 080 
°- ±zo -0.054 

5.65l 22 3 

1 oo+°- 85 


Q m 

Age[Gyr] 
Ho 


o 26^+°-° 2ti 

U.ZUJ„ 032 

12.9 + H 
76.51$ 


0.266™ 
13.211;° 
74.5l|? 


U.ZUO_ 030 

13.111;° 
75.21?;° 


2fiO+°- U28 

U.ZUU_ 033 
12.1111 

82.0™ 



Table 2: Mean values and 68% C.L. errors from current observations of CMB for the cases 
with WMAP5 alone and all data combined. 



increases. In this case, N„ cannot be less than the standard value of 3.046. Thus it may 
be appropriate to study adopting the prior N u > 3.046 to constrain a scenario with such 
an extra radiation component. We denote it as AN^ = N v — 3.046. With this prior, we 
obtain an upper bound on an extra radiation component as N u < 8.19 (or AA^; xt < 5.14) 
at 95 % C.L. when Y p is taken as a free parameter and improve to be N v < 6.35 (or 
AA^; xt < 3.30) at 95 % C.L. when the BBN relation is assumed. Notice that these limits 
are weaker than those with no priors on N u , which are N u < 7.94 (AA^ xt < 4.89) and 
N v < 6.38 (AA"^ xt < 3.33), respectively. This somewhat peculiar fact stems from the shape 
of the likelihood shown in Fig. [2J, which is not symmetric with respect to N u = 3.046. Since 
the differences due to the N u prior are not negligible, caution is needed when we use these 
constraints regarding the prior on A,,. 

Even when we limit ourselves to the case with three active neutrino species, a deviation 
from the standard value of N u = 3.046 is possible. In a scenario with low (MeV scale) 
reheating temperature, N v can be less than 3.046. In this case, N v only takes the value 
less than the standard one. Thus it may be interesting to investigate a constraint on N v 
assuming N v < 3.046. As regards the treatment of Yp, we do not consider the case with 
adopting the BBN relation here because, in a scenario with MeV reheating temperature, Y p 
should be calculated taking into account the non-thermal neutrino distribution functions 
and oscillation effects [3-6]. These effects drive Y p to increase as N v decreases contrary 
to the usual case where N v just represents a measure of the expansion rate. (This is 
why we are not showing constraints for the prior N v < 3.046 with the BBN relation in 
Table [3j) Since taking into account this effect is beyond the scope of this paper, we 
show the constraint for the case with Y p being varied freely, which can be considered 
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No priors on N v 


N v > 3.046 


N v < 3.046 




Mean 


68% T 
68% | 


95%T 
95% | 


68%T 
68% | 


95%T 
95% | 


68%T 
68%1 


95% t 
95%| 


CMB all 

(Y p . free) 


4.24 


5.47 
2.03 


7.94 
0.96 


5.51 
(3.046) 


8.19 
(3.046) 


(3.046) 
2.05 


(3.046) 
1.17 


CMB all 

(Y p : BBN relation) 


3.71 


4.80 
2.27 


6.38 
1.39 


4.70 
(3.046) 


6.35 
(3.046) 




CMB all 

(Y p = 0.24: fixed) 


3.89 


4.89 
2.19 


6.84 
1.28 


4.87 
(3.046) 


6.88 
(3.046) 


(3.046) 
2.12 


(3.046) 
1.27 


Planck 
(Y p . free) 


3.11 


3.44 
2.72 


3.83 
2.41 


3.45 
(3.046) 


3.87 
(3.046) 


(3.046) 
2.72 


(3.046) 
2.43 


Planck 

(Y p : BBN relation) 


3.06 


3.26 
2.87 


3.44 
2.68 


3.25 
(3.046) 


3.44 
(3.046) 




Planck 

(Y p = 0.24: fixed) 


3.19 


3.43 
2.95 


3.67 
2.72 


3.37 
(3.046) 


3.63 
(3.046) 


(3.046) 
2.87 


(3.046) 
2.67 



Table 3: The mean values and 68% and 95% limits of N v for current and future CMB 
data. 



as the conservative one for the prior N v < 3.046. We obtained constraints N v > 1.27 
and N v > 1.17 at 95% C.L. for the cases with Y v being fixed as Y p = 0.24 and Y p 
being assumed as a free parameter, respectively. For a scenario with MeV scale reheating 
temperature, these limits are translated into the lower bound on the reheating temperature 
as T rch > 2.0 MeV [6]. 

Finally, we investigate a future constraint on cosmological parameters paying particular 
attention to N v and its effects on constraints on other parameters. We use the data 
expected from the future Planck experiment and make a MCMC analysis following the 
method in Ref. [57]. As for the specification of Planck, we adopt the following parameters 
for the instrument. For the frequency channels of v — 100, 143 and 217 GHz, the width 
of the beam and the sensitivities per pixel for temperature and polarization are adopted 
as (^ FW Hw[arcmin],a T [/iK],cTp[/iK]) = (9.5, 6.8, 10.9), (7.1, 6.0, 11.4) and (5.0,13.1,26.7), 
respectively. Other frequency channels are assumed to be used to remove foregrounds. 
We make use of the data up to I — 2500 in order that our results will not be affected by 
the SZ effect and the marginalization over A$z is not performed. These setups for the 
Planck forecast are similar to the recent works performed in Refs. [26,30,51], but explored 
parameter spaces are different. We make a simple extension by adding N v and Y p to 
the standard 6 dimensional parameter space, but theirs include neutrino masses and/or 
lepton asymmetry. When one would like to check the constraint on an extra radiation 
component in a simple scenario, one can refer our results here. However, when some 
other particular setups are considered such as a scenario with large lepton asymmetry and 



11 



massive neutrinos, the above mentioned references should be consulted. 

Our results are summarized in Tables [3] and HI As seen from the Tabled the constraint 
is most stringent when the BBN relation is adopted, and in this case, we obtained a future 
constraint as 2.68 < N v < 3.44 at 95 % C.L. Another point which should be noted is that 
fixing of Y p — 0.24 can bias the determination of some other cosmological parameters such 
as Uf, and n s , which was already pointed out in Refs. [51,52]. However, when we vary the 
value of N u , the effect of fixing of Y p = 0.24 is partly cancelled by the change in N v . In 
fact, this in turn results in biases of other cosmological parameters such as u c and 9 S which 
are strongly correlated with A^. Therefore, Y p should be carefully treated in investigating 
cosmological constraints with future CMB data. 



parameters 


Planck 
(Y p . free) 


Planck 
{Y p : BBN relation) 


Planck 
{Yp = 0.24) 


uj b 

0s 
T 

n s 

lntlO 10 ^) 

Y P 
N v 


u.uzz,io_ 0002g 

n 1 1 ns+ 00046 

U.llU8_ 00056 

1 n4n4+oooi4 

1.U^U^_q 0014 

n nssi +0.0050 
0R4+0.009 

o 066 +a016 
o.uoo_ 016 

246 +a020 

q i i +0.33 
°- ±i -0.39 


n n997c;+ u -ouo26 

U.UZZ / U_ oo27 

noi +a0028 

U.llUl_o 0028 

1 04060+ 000044 

±.U^±UUU_o .00049 

n nssi+o- 0053 

U.USSl_o.o063 

Qfi4+o o 10 

u ' yu ^-0.010 

3-065±g : gJJ 
0.248818X 
3 06 +0 - 20 

°- uu -0.19 


n fi997^+ u -00027 
U.U^^(O_ 00026 

1 1 90+ ' 0033 
U.±±ZU_ 0036 

i.o4ooo^;Zl 

0.0880l°S 
0.963±g;gg 

3.068jS;gg 

1 q+0.24 
°- ±J -0.24 


Age[Gyr] 
H 


0.256±8Z 
13.63185* 
72.3+|4 


o.256±8:gss 

13.67^; 20 
72-OlhS 


n oc;c;+0.009 
u.zoo_ 010 

13.561*1 
72.7111 



Table 4: Forecasts on mean values and 68% errors of A^^ and other cosmological parameters. 

We would like in the end to comment on how our discussion so far can be affected 
by theoretical uncertainties in the recombination process [58-66] . Since the change of Y p 
can influence the recombination process, its uncertainties might affect the cosmological 
parameter determination in some way. Thus it may be worth mentioning here on the 
effects. For this purpose, we proceed with the same analysis as have been done in Ref. [52] 
but varying N u here. Two parameters Fjj and fe# e , which represent the uncertainties in 
the recombination modeling, are included among other free parameters. (See Ref. [52] 
and references therein for more details). We impose top-hat priors, < Fh < 2 and 
< hue < 1-5, which are very conservative ones, to take into account the uncertainties 
in the recombination theory. We made the analyses for the two cases where Y p is given 
from the BBN relation and Y p is treated as a free parameter. In both cases, we found that 
the constraints on other cosmological parameters including N v are scarcely affected even 
by very conservative prior on Fh and bn e - The mean values are unchanged and errors 
increase only very slightly (no more than 10% for any parameters other than Fh and bn e )- 
Therefore we can say that the uncertainties parametrized with F H and bn e do not change 
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much our results of the Planck forecast discussed above. However, we would need more 
understanding of uncertainties in the recombination theory for the precise determination 
of cosmological parameters in future CMB surveys. 

4 Summary 

We discussed the issue of probing the effective number of neutrino species N v from CMB 
data alone. Although a constraint on N v has been investigated by many authors, in most 
analyses, some combinations of data set such as CMB+LSS, CMB+ifo; CMB+LSS+ifo 
have been used to constrain N u . This is partly because N u has severe degeneracies in 
WMAP with some other cosmological parameters such as u m and h; thus N v can be more 
constrained by combining some data sets. However, when we combine data from LSS, 
some subtleties can arise: a constraint from LSS data depends on how we treat non-linear 
corrections/bias. Furthermore, different galaxy data lead to slightly different constraints 
on N u . In addition, as for the Hubble prior, the prior usually adopted is H = 72 ±8 which 
is from the result of Freedman et al [31]. However, another group has reported somewhat 
different value as H = 62.3 ± 5.2 [32]. Since different priors on the Hubble constant can 
give different results, in this respect, the constraint obtained by assuming some prior on 
Ho should be regarded taking into account the above uncertainty. Taking these issues into 
consideration, it may be interesting to study a constraint on N v removing such subtleties, 
which can be done by using CMB data alone. 

In this paper, first we discussed the effects of N u on CMB and the issues of degeneracies 
with some other cosmological parameters. Phenomenological descriptions of its effects 
on the heights and the positions of acoustic peaks were also given. Then, in section [31 
a constraint on N v was studied by using CMB data alone. We made use of the data 
from WMAP5, ACBAR, BOOMERANG and CBI. As discussed there, although WMAP 
measurement is very accurate, its precision is limited up to the 2nd peak/3rd peak. We 
have explicitly shown that the information up to the 2nd/3rd peak is not enough to 
constrain N u severely. This was demonstrated by making the analysis with WMAP data 
alone, in which a sensible constraint cannot be obtained. However, if we include the data 
on small scales, the degeneracies of N v with some other cosmological parameters can be 
removed to some extent; then a stronger constraint can be obtained. In fact, on small 
scales, the amplitude is suppressed due to the free streaming effect by increasing N v , which 
is similar to the effects of Y p through the diffusion damping. Thus we have studied the 
constraint on N v assuming different priors on Y p : adopting the BBN relation to derive 
Y p for given N v and u b) assuming 7 p as a free parameter and usual fixing of Y p = 0.24. 
Depending on the prior, the constraint slightly changes. We obtained the 95 % limits as 
0.96 < N u < 7.94 for the case with Y p being free, 1.28 < N u < 6.84 for Y p being fixed 
as Y p = 0.24 and 1.39 < N v < 6.38 when the BBN relation being adopted. It should be 
noted that these constraints are comparable to that obtained using CMB+LSS in previous 
works. 
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One of the main purposes of constraining the effective number of neutrino species using 
cosmological data is to check the standard value of N v independently from particle physics 
experiments. Thus we primarily focus on the analysis with no prior on N u . However, from 
the viewpoint of constraining extra radiation which may be motivated from some particle 
physics models, a constraint obtained by assuming the prior N v > 3.046 may be interesting 
since an extra radiation always increases the value of N u . In this respect, we also made 
an analysis adopting this prior and obtained the constraint on the effective number of 
neutrino as N v < 6.35 and N v < 8.14 for the cases where the BBN relation is adopted and 
Y p is treated as a free parameter. 

On the other hand, in a scenario with low-reheating temperature, the effective number 
of neutrino species can be reduced. In this case, another prior may be motivated to be 
assumed for a simple scenario of low-reheating temperature with three relativistic neutrino 
species. In this regard, we have also studied the case with the prior N v < 3.046 and 
obtained the constraints as N v > 1.17 for the cases with Y p being assumed as a free 
parameter. This can be translated into the lower bound on the reheating temperature as 
T reh > 2.0 MeV. 

We have also discussed a future constraint on N v using the expected data from Planck 
experiment. It was shown that the attainable constraint on N v from Planck is 2.68 < 
N v < 3.44 at 95% C.L. when the BBN relation is adopted for Y p , which is most stringent 
compared to the other cases. Since Planck experiment can probe CMB down to smaller 
scales than WMAP, Planck alone can give a stringent constraint on N v . 

The interplay between particle physics and cosmology is now becoming more important 
in the era of precision cosmology. One of such examples is the number of neutrino species, 
which was investigated in this paper. In light of upcoming more precise observations 
of cosmology, research of this kind will bring us fruitful insight for particle physics and 
cosmology. 
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